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Abstract 

We present here a convergence proof for spectral approximations for 
hyperbolic systems with initial and boundary conditions. We treat in detail 
Chebyshev collocation, but the final result is readily applicable to other 
spectral methods, such as Legendre collocation or tau-methods. 
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INTRODUCTION 


In the paper [1], we derived stability results for spectral methods 
applied to initial-boundary value problems for hyperbolic systems. The paper 
demonstrates that one can bound certain weighted spatial norms of the 
solution in terms of norms of the boundary data (homogenous initial conditions 
are assumed). The bounds also contain powers of N, which is the degree of 
the approximating polynomials. 

Here we show that the approximations discussed above actually converge to 
the exact solution, at least when this solution is smooth. We bound the error 
in the numerical method by a power of N multiplying a term which depends 
only on the exact solution — more precisely, this is the interpolation error of 
the initial value and boundary derivatives. For sufficiently differentiable 
functions, this interpolation error will decay fast enough to drive the full 
approximation error to zero. We have not attempted to derive the sharpest 
bound of this type, but merely to show that such a bound exists. 

The method of proof here is similar to the one in [1,2], where basic 
results are first deduced for a scalar equation, and then extended to the full 
system. Accordingly, the paper is divided into two sections, the first 
dealing with the scalar case and the second with the system. By the means of 
Gauss-Lobatto quadrature formulas, we first bound the error at outflow for a 
single scalar equation. Then we use this estimate, together with the basic 
stability result of [1], to bound the overall error for a system. 
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1. THE SCALAR CASE - CHEBYSHEV COLLOCATION 

Given the equation 


( 1 . 1 ) 


! u t = u x |x| < 1, t > 0 

u(x,0) = f(x) 
u(l,t) = g(t). 


we consider the pseudospectral method based on collocating at the extrema of 
T N+1» w h ere T m is the Chebyshev polynomial of degree m 


T (x) = cos(m cos *(x)) . 


It is shown in [3,4] that the pseudospectral approximation v = v^(x,t) 
satisfies 


/ 


9v 9v 
9 1 9"x 


x(t) 


T' 

N+l 


(x) 


( 1 . 2 ) 


v(x,0) = 


P N f 


v( 1 , t) = g(t). 


with the projection operator Pfj defined by 


f is a polynomial of degree < N 

(1.3) 

( P N f )( x ) = f ( x ) at x = x j = cos , 


0 < j < N 
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Define now 6^(x,t) as the difference between vjj and the projection of u 


<S N (x,t) = v - P N u. 


This function satisfies the equation: 


(1.4) 


< 3T S N <X ’ t) 


\ « N (x,0) 


37 V *’ 0 


f P N 3x 3x P N + T T N+1 


= 0 


V « N (i»t) = °* 


We note at this stage that the polynomial (l+x)6 N T' + ^ is of degree 2N+1 
and therefore may be integrated exactly by the Gauss-Lobatto quadrature 
rule. This results in the following 


.1 T ' +1 


N+l 


(1.5) 


"I J i 2 

/ 1 - x 


= ? I H7 (1 + > W W’V ‘ °> 

j=0 j J J J 


C 0 = C N+1 = 2 » c j = 1 for 0 < J 1 N, 


since 1 + x N+1 = 0, 6 N (x Q ) = 0 and Tj^ +1 (xj) = 0 for j = 1,**.,N. In 

fact, because of the term (1-x) and the boundary condition, the indices in 
the quadrature sums may run only from 1 to N (instead of 0 to N+l). 

Defining 

= ITx ^ P N " P N Sx ’ 


multiply (1.4) by (1 + x)6^ and integrate to get 
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1 

2 dt 


1 (1+x)6 n (Xjt) 

/ — — dx 


-1 


/ 1 - 


1 (l+x)S N (x,t) 9 

-1 ~ ' 2 6 

1 / 1 - x 2 


( 1 . 6 ) 


1 (l+x)6 N (x,t)Q N 

+ / — dx. 

-1 / , 2 
/ 1 - x 


An integration by parts on the right-hand side produces: 


1 d 1 ( 1+x ) 5 N ( x >t)dx 

2 dt / — ■ dx ° 


-1 


7 / 


6^(x,t) 


dx 


✓ 1 - 


* (1-x) / 1 - x^ 


(1.7) 


+ J 

-1 


1 (l+ x )<5 N ( x ,t)Q N ( X ,t) 


✓ 1 - 


We again use the Gauss-Lobatto formula and reach 


1 d r 1 ( 1+x ) 6 N ( x »t) 

2dF J — dx 


-1 / . 2 
V 1 “ X 


( 1 . 8 ) 


N 6 (x. ,t) 

7T N j * 

2N / 1 - x. 

3 = 1 J 


8N 6 N^ -1, t) +j|- l (1 + x )6 N (x. 

3 = 1 J J 


< - 


N 6 (x.,t) 
tt v N j 9 ' 

2 n ^ nr - x. 

3 = 1 3 


_ o N 6 (x.,t) 

TT v 7T ^ N_J ’ 


— 6 f-i t ) + — y — J 

8N °N 1 l ' t} + 2N .<■ 1 - x. 

3 = 1 3 


N 


+ I (1 + X P (1 " x.)Q M (x,,t). 


3 = 1 


3 N j 


(x,t)dx 


dx. 


.t)Q N (Xj ,t) 
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Equation (1.8) yields immediately — since 6,,(x,0) =0 — 


(1.9) J 


1 (1+X) . 2 / ,2 


-1 


/l-x 


6 2 (x,t)dx+^ / ' 6jJ(-l,t)dt < J |||Q(x,t)||| dt 
2 w ™ 0 N 0 


where 


( 1 . 10 ) 


| Q(x, t) | | | 2 -i ^ (1+x.)^ (l-x.)Q^(x.,t: 


In the next section we will need a different version of (1.9), in which 
the time integral is weighted by e , n > 0: 


(1.11) / e 2nt / 6 2 (x,t)dxdt + / e 2nt 6 2 (-l,t)dt 

0 -1 " 2 N n 0 


/l- 


x 


< 2^ / e 2nt I I lQ N C X ,0| | I 2 dt. 

This form matches the Laplace—Fourier transforms which are used in the basic 
stability estimates. 

o 

It should be noted that the bound on 5 N (-l,t) obtained in (1.11) 

(1.12) J e" 2nt 6 2 (-l,t)dt < cx(N) / e" 2nt | | |Q N (x,t)| | | 2 dt. 


4N 

a(N) = H 


is very crude. A better estimate can be found for the boundary error u(-l,t) 
- v(-l,t) by taking the Laplace transform of (1.1) and (1.2) and analyzing 
the difference. This approach was used by Dubiner [5] who obtained 
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(u v)( l»t) ~ N . The same bound holds for P N u - v, since, for smooth 
functions, P N u approaches u faster than any power of i . 


2, CONVERGENCE PROOF FOR SYSTEMS 

Consider the first-order hyperbolic system of partial differential 
equations 


( 2 . 1 ) 


9 u _ . 9u 

at A ax ’ 


x| < 1, t > 0 


where 

u = u(x,t) = (u (1) (x,t), u (2) (x,t),..., u (n) (x,t))' 

is the vector of unknowns and A is a fixed nx n coefficient matrix. Since by 

hyperbolicity A is similar to a real diagonal matrix we may, without loss of 
generality, take it diagonal: 



The solution of this system is uniquely determined if we specify initial 


conditions 
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(2.3) 


u(x,0) = f(x) = (f I (x), f rI (x)) 


and boundary conditions 


(2.4) 


i uVl.t) = Lu II (— 1 , t) + g^t) 
u n (l,t) = Ru I (l,t) + g U (t). 


In these formulas, f and 
vectors, and 


(2.5) 


g = g(t) = (g I (t), g IX (t)) are prescribed n- 
( (»\ u 11 ) 

{(f 1 . f 11 ) 

g 11 ) 


is the partition of these vectors into inflow and outflow components — 
corresponding to the partition of A in (2* lb). L and R are constant 

reflection matrices of order l x (n-£) and (n-Jl) x respectively. 

We shall discuss only problems whose solutions decay in time and 
therefore postulate: 

ASSUMPTION I. There exists a constant y > 0, such that 


|r| |l| < l - y < l 
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( I A | and |v| denote the Euclidean norm of a matrix A or vector v, 
respectively.) 

In a pseudospectral Chebyshev approximation to (2.1), one seeks a vector 
v = vjj of N-degree polynomials such that 


(2.6a) 

The multiplier t 


(2.6.b) 


3t A 3x + T N+1^‘ 


T^) is determined by the boundary conditions 


v (-l,t) = Lv II (-l , t) + g T (t) 


» n (l,t) = Rv I (l,t) + g I1L (t). 


Initially, v is defined by collocation: 


(2.6c) 


v I (x,0) = f I (x) 
v II (x,0) = f II (x) 


at x = x.. , j = 1 ,2, • • • ,N+1 

at x = Xj , j = 0,1,* •• ,N. 


In [1] the stability of the approximation (2.6) has been established, under 
zero initial conditions. Here we shall prove the convergence of v^(x,t) 
to u(x,t) as N tends to infinity. 

We define a pair of projection operators P= (p 1 , p 11 ) by requiring that 
for any function F, P* F and P^ F be polynomials of degree N at most. 


satisfying: 



(2.7) 


(P F)(Xj ) = F(Xj ), 


(P 11 F)( Xj ) = F( Xj ), 


j = 1 »2,« • • ,N+1 


j = 0,1,... ,N. 


Note that two distinct sets of nodes are used for collocation - cf. (2.6c). 
We can state now the main convergence result: 


THEOREM. Let e = e^(x,t) = v-Pu be the error in the pseudospectral 
approximation (2.6) to the hyperbolic system (1.1). Let Q = (Q 1 , Q 11 ) be 
the approximation error 


( 2 . 8 ) 


Then 


i 1 - ‘V 1 lr - k pI - 1 ) 


n . . |_ P IX n } 

^ 3x 3x J 


(2.9) 


/ e 

0 


" 2nt " e ( x , t ) II 2 dt < - K - N - ” (N) / e~ 2nt | | | Q(x, t) | j | 2 dt 


KNcx(N) 


2n 


with K independent of N, a(N) as defined in (1.12) and 
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( 2 . 10 ) 


lle(x,t)ll 2 = / (1-x) |e I (x,t) | 2 — . dx . + / (l+x) |e I]: (x,t)| 2 

-1 , . 2-1 

Y 1-X 

| | |Q(x,t)| | I =1 l (1-x ) 2 (1+X ) |Q I (x. , t ) | 2 
j=l J J J 



N 2 TT . 

I (l+x.) (1-x ) |Q (x ,t)| . 

j-1 J J J 


(Here we generalize the seminorm defined in (1.11).) 


Proof: Let r be the solution of 


( 2 . 11 ) 


3r _ 8r 

at ax 

r(x,0) = u(x,0) 
r I (-l,t) = u I (-l,t) 
r n (l,t) = 


and let s be the pseudospectral Chebyshev approximation to (2.11), i 


3 s .3s _ 

at “ A ax + t n+i^ 

s(x,0) = Pu(x,0) 
s I (-l,t) = u I (-l,t) 
s I]C (l,t) = u II (l,t). 


dx 



• 6 • * 


( 2 . 12 ) 
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lt is obvious that r = u, but s does not satisfy the boundary conditions 
(2.6c) and therefore its multiplier 0 is distinct from t. In any case, for 
6 = s - Pr, we have: 

(2.13) 

/ e“ 2nt 115 (x,t) II 2 dt + / e“ 2nt |5 I (1 ,t) | 2 dt 

0 n 0 

+ m f e_2nt l 6lI (-l.t>|dt < e“ 2nt | | | Q(x, t) | | | 2 dt . 

^ n 0 ^ 0 

This is clearly a restatement of (1.11). Next, we compare the spectral 
solution s with v — as defined by (2.6). We can show that 

(2.14) / e" 2r,t II s-vll 2 dt < - N J e _2r,t [ |6 I (1 ,t) | 2 + | 5 XI (-1 , t ) | 2 ]dt . 

0 11 0 

Indeed, v-s satisfies: 

ft <v_s) * k <v - s) + T m-i (, ' 6) 

(v-s)(x,0) = 0 

(v 1 - s I )(-l,t) = L(v n - s II )(-i;t) + Ls II (-l,t) - u I (-l,t) + g 1 
= L(v TI - s II )(-l,t) + L6 IX (-1 ,t) , 

(v 11 - s II )(-l,t) = R(v* - S 1 )(l,t) + R6 I (l,t), 
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and the inequality (2*14) follows from the stability estimate in [1, Theorem 
5.1]. 

Now, since Pr = Pu we have 


e - 6 = v-s, 

and, hence, 

00 00 00 
/ e~ 2nt lie (x,t)l! 2 dt < J e _2nt Il6(x,t)0 2 dt + / e' 2nt II v-s II 2 dt. 

0 0 0 

The first term is majorized in (2.13) and the second in (2.14), thus 
establishing the theorem. 

Again, we emphasize that the bound for u — v^j is the same as the bound 
for Pu - Vjj which we just have computed — for smooth functions u, which 
Pu approaches rapidly. 

We conclude with two remarks: 

a. The stability estimate of (2.14) explicitly uses boundary values. 
This is why these have to be bounded beforehand, which was done in the 
preceding section, formula (1.10). 

b. Our result applies not only to Chebyshev collocation, but also to 

other spectral methods (which, however, have to satisfy assumptions I, II and 

III of [1]). Indeed, the only quantity that varies with the spectral method 

employed is the coefficient weight at ±1 for the Gauss-Lobatto quadrature. 

Once it is known that this weight is ~ -L- —which is the case for Gegenbauer 

N 

collocation, as shown in [1]— the same proof follows through. In particular, 
it is sufficient to evaluate outflow errors for a scalar equation in order to 
estimate errors at both boundary points for a system. 
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